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SUMM/LRI 


This  paper  presents  a  proposed  set  of  acceptance-sainpling  plans  for 
life  testing  and  reliability  when  the  underlying  life  distribution  is  of 
the  WeibuU  form*  Inspection  of  the  sample  is  by  attributes  with  the 
life  test  trvincated  at  a  preassigned  time^  t«  A  set  of  conversion  tables 
is  also  provided  from  which  attribute  sampling-inspection  plans  of  any 
desired  form  may  be  designed  for  the  WeibuU  model  or  from  which  the 
operating  characteristics  of  any  given  plan  may  be  determined.  A  pro¬ 
cedure  using  these  tables  for  applying  the  MIL-STD-105B  plans  to  reliabili¬ 
ty  and  life-testing  applications  is  included. 
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INTRODUCTION 


The  paper  is  a  generalization  of  papers  by  Sobel  and  Tischendorf*- 
and  by  Epstein  that  appeared  respectively  in  the  Proceedings  of  the 
Fifth  and  Sixth  National  Symposium  on  Reliability  and  Quality  Control  in 
Electronics.  Related  work  has  also  been  done  by  Gupta  and  Groll^  who 
have  extended  the  Sobel  and  Tischendorf  procedures  from  the  exponential 
form  to  the  gamma  form.  The  gamma  variable  is  the  sum  of  exponential 
variables  and  hence  the  exponential  model  is  a  special  case  of  the  gamma. 
In  the  two  papers  first  cited  the  authors  assume  that  the  underlying  life 
density  is  exponential  (Eq.  1)  whereas  in  this  paper  the  WeibuU  form  is 
assumed  (Eq.  2).  The  exponential  distribution  is  a  special  case  of  the 
Weibull  and  so  will  be  covered  in  the  plans  and  conversion  tables. 

f(x)  =  (l/n)  exp  [-x/^1,  p  >  0^  X  >  0  ^1  j 

f(x)  =  (P/t)  )(x/r]  f'  '  exp  [-(x/q)^  q  >  0,  x  >  0  (2) 

Both  f(x)  are  equal  to  zero,  otherwise.  In  these  and  the  equations  that 
follow,  X  is  a  random  variable  which  represents  item  life  for  which  x 
is  its  value,  [i,  represents  mean  item  life  for  the  population,  and  p 
is  the  symbol  for  the  shape  parameter  for  the  Weibull  distribution.  For 
simplification  in  discussion  and  computation,  the  characteristic  life, 
n  has  been  used.  For  the  Weibull  distribution  ' 

n  =  p/r(l/p  +  1) 

For  further  discussion  of  the  Weibull  distribution  as  a  statistical  model 
for  lifelength  of  components  or  systems,  reference  may  be  made  to  a  paper 
by  Kao^  in  the  Proceedings  of  the  Sixth  National  Symposium  on  Reliability 
and  Quality  Control  in  Electronics. 


From  Eq.  (3)>  it  will  be  noted  that  for  a  given  q  the  WeibuU  nean 
life  depends  on  the  shape  parameter,  p .  The  effect  of  differences  in  the 
parameter  on  the  shape  of  the  distribution,  as  well  as  the  general  nature 
of  the  WeibuU  model,  may  be  observed  study  of  Figure  1.  In  this 
figure  the  WeibuU  probability  density  function  has  been  plotted  for 
various  value  of  p.  A  plot  for  P  ••  1,  the  exponential  case,  has  been  in¬ 
cluded  for  reference.  This  initial  set  of  WeibuU  sampling  plans  is  for 
product  for  which  the  value  of  this  parameter  is  known  or  can  be  assumed 
to  approximate  some  given  value.  Conversion  tables  and  sampling  plans 
are  provided  for  nine  values  for  p  ranging  from  1/3  to 

A  relatively  small  number  but  a  broad  range  of  values  for  P  was 
selected  for  this  Initial  study.  The  principal  objectives  were  to  develop 
practical  methods  and  techniques  and  to  explore  the  effect  of  differences 
in  value  for  this  parameter  (which  is  the  key  one  for  the  WeibuU  distri¬ 
bution).  As  the  use  of  this  distribution  as  a  statistical  model  Increases, 
additional  conversion  tables  and  sampling  inspection  plans  may  be  con¬ 
structed  for  intervening  values  for  P,  particularly  in  the  widely  encoun¬ 
tered  region  ranging  from  l/2  to  2. 

For  the  procedures  and  plans  developed  in  this  study,  lot  or  product 
quality  is  evaluated  in  terms  of  mean  item  life,  p..  Subsequent  work  has 
been  planned  in  which  related  conversion  tables  and  sampling  inspection 
procedures  will  be  developed  for  application  when  lot  quality  must  be 
evaluated  in  terms  of  the  instantaneous  failure  rate,  Z(t),  at  some 
specified  life  or  future  time,  t. 

FORil  OF  ACCEPTANCE  CRITERIA 

For  the  plans  considered  in  this  paper,  the  foUowing  acceptance¬ 
sampling  procedure  for  life  testing  has  been  assumed: 
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1.  Select  a  random  sample  of  n  items  from  the  lot. 

2.  Put  the  sample  items  to  life  test  for  some  preassigned  period 
of  t  time  units. 

3«  Denote  by  y  the  number  of  failures  observed  prior  to  time  t. 

li.  Accept  the  lot  if  y  S  c,  the  acceptance  number;  if  y  >  c, 
reject  the  lot. 

Curtailed  Inspection  for  accepting  lots  prior  to  t  is  possible  for  the 
rejection  of  the  lot  since  it  is  possible  to  observe  (c  +  1)  failures 
before  time  t. 

Note  that  this  acceptance  procedure  is  the  same  as  that  specified  for 
the  MIL-STD-105B^  san^Jling  plans  with  the  exception  of  the  introduction 
of  a  testing  truncation  time,  t.  It  is  also  possible  (as  for  the  105-B 
plans)  to  employ  double  or  multiple  sampling  instead  of  single  sampling 
as  described  above  and  by  so  doing  reduce  the  average  number  of  items 
at  pt  -  AQL  that  must  be  put  on  life  test.  However,  the  "economy"  achieved 
is  at  the  expense  of  longer  elapsed  testing  time. 

The  probability  of  acceptance  for  a  lot,  P(A),  under  plans  of  the 
above  form  depends  on  the  probability,  p',  of  item  life  being  less  than 
(or  equal  to)  the  test  truncation  time,  t.  For  cases  for  which  p  is 
known  and  with  time,  t,  preassigned,  p«  is  thus  a  function  of  mean  item 
life,  p,  only.  The  operating  characteristics  of  any  specified  plan  thus 
depend  only  on  t  and  p.  In  order  to  provide  tables  for  general  use  in 
the  design  or  evaluation  of  plans  for  any  application  rather  than  working 
in  terras  of  specific  values  for  t  and  p,  the  dimensionless  ratio  t/p  will 
be  used.  In  application  of  the  plans  or  tables  to  a  specific  application, 
conversion  between  the  ratio  and  specific  t  and  p  values  is  extremely 
easy. 
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A  Of  tables  has  been  computed  to  provide  for  the  Wei. 

bull  distribution  the  connection  between  the  dimensionless  quantity  t/p 

P  (Tables  1  and  2).  with  these  tables  acoeptance.saiiipling  plans  of 

desired  fo:m  can  be  designed  or  evaluated  using  attribute  sampling  theories 
and  practice. 

For  cases  for  which  the  lot  site,  w.  i.  large  in  relation  to  the 
sample  size.  n.  the  number  of  failures  prior  to  t  approaimtes  the  bi- 
uomial  distribution  with  para,«ters  n  and  p.,  where  p.  is  defined  as  the 
area  under  the  life-length  distribution  up  t„  t.  The  probabUity  of 

acceptance  P(a,  depends  on  the  cumulative  number  of  faUu.es  prior  to  time 
t.  This  probability  is  given  by 

P(A)  =  P  (y  S  c)  =  2  p.y  (,.p.)n-y 

y=o  (^) 

The  binomial  distribution  has  been  assumed  for  the  sampling  plans  given  in 

this  paper  e»=ept  for  cases  for  which  the  sample  size  is  relatively  large. 

For  these,  the  Polsscn  distribution  has  been  used  as  an  approximation  to 

the  binomial.  The  probabUity  of  acceptance  for  the  Poisson  is  given  by 

P(A)  =  P  {y  S  c}  =  £  -Lnp*)^  g-np» 

y=o  y!  (5) 

An  important  potential  use  for  the  conversion  tables  provided  in 
this  paper  is  in  the  adaptation  of  the  MIL.SIC-105B  plans  to  reliability 
and  Ufe-testing  applications.  In  describing  the  operating  characteris- 
tics  of  these  IOSb  plans,  the  quality  of  submitted  lots  is  measured  in 
terms  of  pi,  the  per  cent  defective.  With  the  conversion  factors  this 
form  Of  description  may  be  converted  directly  to  measurement  in  terms  of 
the  t/p  ratio.  With  this  conversion  the  lOSB  plans  may  be  cataloged  for 
appropriate  choice  in  reliability  applications.  (Plans  have  been  made 
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for  the  preparetion  of  tables  listing  for  various  values  of  e  this  ratio 
at  the  AQl  and  the  MPD  for  each  plan  in  the  105B  nanual.)  Alternatively, 
if  some  lOSB  plan  has  been  selected,  its  operating  characteristic  curve 
-ay  be  determined  in  terms  of  the  t/p  ratio,  or  if  the  testing  time,  t, 
been  specified,  in  terms  of  the  lot  mean,  p.  An  example  eng>loying 
such  a  conversion  is  shoun  later  in  this  paper.  It  should  also  be  noted 
that  with  the  matching  plans  provided  in  the  lOSB  collection,  the  options 
of  double-sampling  and  multiple-sampling  are  also  available.  The  --pn- 
sises  and  acceptance  numbers  listed  may  be  used  and  the  established  pro¬ 
cedures  for  employing  this  form  of  sampling  i„  attribute  Inspection  may  be 
followeda 


COMPUTATION  OF  CONVERSION  TABLES 


The  probability,  p«,  of  an  item  failing  before  the  end  of  test  time 
t  IS  given  by  the  cumulative  distribution  function  (c.d.f.).  For  the 
Weibull  model  the  equation  for  this  function  is 

c.d.f.  o  F(x)  ■  1  -  ejqj  [-x^/a]. 

The  equation  for  the  mean,  n,  of  the  Weibull  distribution  is 

M.  -  r(i/p  +  I)  . 

Computations  may  be  simplified  by  the  following  substitutions: 


b  =  l/P 


Equation  (6)  then  becomes  ~ 


(8) 

(9) 


F(x)  •  1  -  ejqj  , 


(10) 
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and  Equation  (7)  becomes 


-  n  r  (b+l) 


(II) 


The  probability,  p»,  of  an  item  failing  before  the  end  of  test  time, 
t,  is  thus  given  by 


F(t)  -  1  -  e:qp  (-(t/q) =  p'  . 

This  may  be  rewritten  as 

1/(1  -  P*)  -  exp  ((t/q)'^*^]  , 

which  in  turn  may  be  converted  to 

-ln(l-pi)  -  (t/q)*''*' 

where  In  denotes  the  Naperian  logarithm. 

Raising  both  sides  of  the  above  equation  to  the  b  power  gives 
(“InCl-p')]*^  =  t/q  . 

This  equation  may  be  solved  for  q  to  give 

q  =  t/[-|n(l-p')]^  . 

But  from  Equation  (11), 

q  =  M/r(b+( )  . 

Substitution  of  this  value  for  q  in  Equation  (16)  gives; 
p/r(b+|)  =  t/  [-ln(l-p')]‘’  , 

or  t/p.  »  (-ln(l-p')]^  /  r(b+l )  , 


(12) 

(13) 

(11^) 

(15) 

(16) 

(17) 

(18) 
(19) 


This  equation  establishes  the  relationship  between  the  dimensionless  ratio 
t/p  and  p',  the  probability  of  item  life  being  equal  to  or  less  than  t. 
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It  may  be  noted  that  for  the  attribute  form  of  sampling  inspection 
considered  here,  only  this  dimensionless  ratio  between  test  time,  t,  and 
Item  mean  life,  p,  need  be  of  concern.  The  Weibull  scale  parameter, a 
(or  its  equivalent  characteristic  life,  „  )  has  been  eliminated.  In  the 
mathematics  of  these  plans  and  procedures  it  has  been  assumed  that  the 
Weibull  location  parameter,  y  ,has  a  value  of  0.  If  in  application, 
however,  7  has  some  non-zero  value,  all  that  is  necessaiy  is  to  subtract 
the  value  for  7  from  the  value  for  t  to  get  t^,  and  from  the  true  lot 
mean  p,  to  get  p^.  These  converted  values,  t^  and  p^,  are  then  used  for 
all  computations.  Any  solutions  in  terms  of  t^  or  p^  can  be  readily  con¬ 
verted  back  to  real  or  absolute  values  by  simply  adding  the  value  for  7. 
This  procedure  for  handling  the  location  parameter  will  be  illustrated 

later  in  Example  3.  Only  the  parameter p  (or  b,  which  is  1/p)  must  be 
known. 

To  put  this  relationship  equation  (Eq.  I9)  in  a  form  for  which  nu¬ 
merical  values  for  relationships  may  be  more  easily  computed,  the  fol¬ 
lowing  change  may  be  made: 

t/p  -  f-ln(|-p-)]Vr(b+l)  =  exp  {  In  [- In  (l-p  •  )]  Vr(b+I  ) 

-  exp{  b  In  [-ln(l-p')]J  /r(b+l )  .  (20) 

Values  for  the  expression 

In  (-In  (1  -p.)] 

were  obtained  from  a  table  of  the  inverse  of  the  cumulative  probability 
function  of  extremes  prepared  by  the  National  Bureau  of  Standards.^ 

This  table  tabulates  the  function 

y  -  -  in  (.  1„  .  (25, 
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5y  substituting  (1  -  p‘)  for  9^  the  negative  value  of  Expression  21  is 
obtained.  Values  for  e  raised  to  this  power  were  read  from  the  National 
Bureau  of  Standards  tables  of  the  exponential  function.^  Values  for  the 
gamma  functionj  r(b  ♦  1)^  were  obtained  from  a  table  prepared  by  Dwight.^ 

A  table  of  values  for  the  per  cent  truncation,  (t/p)  x  100$,  for 
various  values  of  p«  has  been  prepared.  It  is  presented  as  Table  1. 
Values  for  p*  range  from  .010$  to  80$  with  the  tabulated  values  selected 
in  accordance  with  a  standard  preferred  number  series.  For  convenience 
in  both  tabulation  and  use,  both  the  ratio  t/p  and  p*  are  expressed  as 
percentages  rather  than  decimal  ratios. 

For  detemaining  without  interpolation  the  value  for  p*  when  some 
rounded  value  for  the  (t/p)  x  100  ratio  is  given,  the  relatively  simple 
task  of  preparing  a  table  of  p*  has  been  carried  out.  Referring  to  pre¬ 
vious  equations  it  was  noted  that 


t/p  -  [-ln(l-p')3‘^/r(b+0  . 


Eq.(20) 


Raising  each  side  of  this  equation  to  the  0  power  gives 

a 

(t/p)  =  -ln(l-p')/[  r(b+l)]P  . 

From  this,  an  expression  giving  the  value  for  p»  is  found.  It  is 

P'  -  1  -  exp(  [  r(b+l)  . 


The  table  of  values  for  p»  for  various  values  for  (t/p)  x  100  is 
presented  as  Table  2.  Values  for  (t/p)  x  100  range  from  .010  to  100. 
Again,  the  values  used  for  tabulation  form  a  preferred  number  series. 
With  this  alternate  table  available  together  with  the  basic  original  one 
(Table  1),  a  conversion  may  readily  be  made  either  way— from  (t/p)  x  100 
to  p*  or  from  p'  to  (t/p,  x  100.  Also,  it  will  be  noted  that  the  two 


supplement  each  other  in  that  3  values  giving  a  compressed  range  of 
figures  in  one  table  give  an  expanded  range  in  the  other.  This  allows 
for  somewhat  more  precise  interpolation  in  conversion.  The  two  together 
supply  basic  data  for  the  design  or  evaluation  of  any  life-testing  and 
reliability  sampling  inspection  plan  based  on  the  Weibull  (or  exponen¬ 
tial)  distribution  and  of  an  attribute  form.  For  general  information,  the 
relationship  between  the  (t/p)  x  100  percentage  and  p»  as  given  by  these 
two  tables  has  been  plotted  in  Figure  2  for  each  of  the  various  p  values. 

ESTIMATION  OF  THL  SHAPE  PARAI'jETER 

In  many  applications  the  shape  parameter,  p,  may  be  known  for  the 
product  in  question.  From  past  analysis  of  life  testing  results,  it  may 
be  established  that  some  value  of  known  magnitude  may  be  expected  regularly 
and  so  may  be  used  in  sampling  inspection  procedures.  For  example,  for 
a  certaxn  class  of  electron  tubes  of  receiving  type,  Kao^  has  found  from 
study  of  approximately  2,000  tubes  of  a  variety  of  types  and  applications 
that  a  value  of  1.7  may  be  appropriate.  For  ball  bearings,  Lieblein  and 
Zelen^°  found  a  mean  value  of  1.5l  with  S0%  of  approximately  5,000  bear¬ 
ings  tested  having  p  in  the  interval  1.17  to  l,7i|» 

For  products  for  which  the  value  for  3  is  not  known,  this  parameter 
must  be  estimated  using  failure  data  from  past  inspection  and  research. 

Such  data  may  be  available  in  either  grouped  or  imgrouped  form. 

Ungrouped  Data 

In  this  case  the  failure  data  will  consist  of  the  exact  life  length 
of  each  of  the  r  items  that  fail  out  of  the  n  that  are  tested.  These 
life  values  may  be  listed  in  order  and  designated  by  the  notation 

0  <  <  <  < 

*2  •  • ' "  *r  -  1  “  \  method  of  maximum  likelihood 
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W  then  he  epplied  and  an  eetinete  for  9,  (g),  obtained  by  solving  the 
following  equation: 


•  (25) 


(l/r)  (  f  ,  („.o,  ^  .  {S5) 

r/p  +  S  In  X. 
i  =  l  I 

Grouped  Data 

For  this  case  the  failure  data  will  consist  of  the  numbers  failing, 
f,  during  each  of  a  series  of  k  conveniently  chosen  inspection  time  in¬ 
tervals,  z.  This  ordered  paired  data  may  be  noted  as  z  ,  f  ;  z  r  . 

1  1'’  2*  2* 

^  -  1’  ^ic  wl^sre  <  Zg  and  where  f +  f ^  + 

....  fj,  =  r.  The  maximum  likelihood  estimate  for  the  shape  parameter 

(and  for  the  scale  parameter,  a,  as  well)  are  obtained  by  maximizing  the 
expression 

4  ■<  M  '"J-I  Ifl 

‘  J-Vr  ’  'J  ^  .  (36, 

For  groupod  data  the  method  of  minimized  ohi-squarea  may  be  used.  Eati- 

mates  for  the  t»o  parametera  are  obtained  by  ml^zing  the  efpreasion 
Z7  -z7  -yP 


lo-  L  i  f  L  i-J  {  ; 

j=l  j=|  J 


-e-^,  . 


Graphical  Method 

The  above  methods  for  estimation,  it  will  bo  noted,  are  quite  in¬ 
volved.  For  aoouraoy  and  eoonomy  in  oomputation  a  high-speed  eleotronic 
computer  meat  be  used.  Hoaever,  a  simple  graphical  method  for  estimation 
of  the  Welbull  parameters  has  been  devised.  Estimates  are  obtained  by 
plotting  failure  data  on  Weibull  probability  paper.  The  method  depends 
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on  the  fact  that  the  Melbull  o.d.f.,  F(x),  given  ly  Eq.  (6)  beconee 

straight-line  equation  upon  a  double  logerithmio  transforation. 
Thus 


r^T(x)  ‘  -In  O  +  p  Jn  Cx)  .  (28; 

This  Weibull  paper  has  In  In  versus  In  coordinates  so  the  c.d.f.  will 
plot  as  a  straight  line.  Convenient  scales  are  provided  for  direct 
plotting  Of  raw  data  and  for  obtaining  the  desired  paraneter  estimates. 

Further  discussion  of  the  above  ostination  methods  may  be  found  in 
papers  by  Kao.'*-  ^  aleo,  in  a  recent  study  by  Weiss  a  method  has  been 

detemined  that  may  be  used  to  estimate  this  parameter  by  transformed 
sample  spacings,^^ 


USE  OF  THE  CONVERSION  TABLES 

One  form  of  application  that  should  be  of  considerable  use  is  that 
of  evaluating  the  quality  protection  afforded  by  a  proposed  or  existing 
acceptance-sampling  plan.  A  possibility  of  immediate  interest  is  the 
use  of  a  plan  from  the  MIL-STD-IO^B  Tables. 

Example  (l) 

Suppose,  for  example,  that  a  lOSB  plan  with  an  AooepUble  Quality 
bevel  (AQL)  of  2.5*  and  Sample  Site  Letter  J  has  been  proposed  for 
use.  Reference  to  the  105B  Tables  shows  that  for  single  e-p-.-g  a 
sample  site  of  75  items  and  an  acceptance  number  of  1,  is  specified. 
Suppose  life  testing  time  is  to  be  80  hours  with  simply  a  count  made 
of  the  test  items  failing  by  the  end  of  that  period.  From  inspection 
experience  with  the  product  to  which  the  plan  is  to  be  applied,  it 
seems  most  appropriate  to  assume  a  Weibull  distribution  with  a  value 
for  B  of  1  2/3.  The  lot  sice  will  be  relatively  large  compared  to 
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the  sample  size  of  75  so  binomial  probabilities  for  sample  items  can 
be  assumed.  Actually,  Table  III  of  MIL-STD-IO^B  specifies  that  the  lot 
size  should  be  from  13OO  to  3200,  $01  to  8OO,  and  I8I  to  300  for 
Inspection  Levels  I,  II,  and  III  respectively. 

The  first  step  is  to  detemine  the  probability  of  acceptance, 

P(A),  for  various  values  of  p».  These  probabilities  can  readily  be 
obtained  from  any  one  of  the  readily  available  tables  of  the  cumulative 
binomial  terms  or  tables  of  the  incomplete  beta  distribution.  They  may 
also  be  read  from  the  operating  characteristic  curves  published  as  a 
part  of  the  MIL-STD-105B  Tables.  A  few  of  these  values  for  this  plan 
are  shown  in  the  first  end  second  columns  of  the  tabulation  below. 

Next,  the  first  of  the  conversion  tables.  Table  1,  is  used  to  obtain 
values  for  the  ratio  (t/n)  x  100  for  each  of  the  p»  values.  These 
table  values  are  listed  in  the  third  column.  Finally,  using  the  value 
for  t  of  80  hours  each  of  the  (t/n)  x  100  ratios  are  converted  to 
values  for  |j..  For  example,  the  ratio  for  a  p«  of  $%  is  18. 8U.  Thus 
(80/p.)  X  100  «>  I808U  or  p  -  hours.  These  computations  have  been 
made  with  results  as  shown  in  the  last  column.  One  may  now  note  that 
if  a  lot  is  submitted  to  this  plan  whose  mean  life  is  215  hours,  the 
probability  of  its  acceptance  is  only  .01  or  one  in  a  hundred^  on  the 
other  hand,  if  the  mean  life  for  a  lot  is  7U5  hours  the  probability 
of  acceptance  is  .98.  These  probability  and  mean  life  figures  based 
upon  t  -  80  hours  can  be  plotted,  if  desired,  to  give  the  operating 
characteristic  curve.  (Of  course  similar  OC  curves  for  other  known 
values  of  t  may  be  plotted.)  This  curve  is  the  one  shown  for  B  ■  1  2/3 
in  Figure  3. 
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'  (in  %) 

P(A) 

(t/n)  X  100 

2 

.98 

10.77 

745 

3 

.92 

13.78 

580 

h 

.82 

16.42 

490 

5 

.68 

18,84 

425 

6  1/2 

•hB 

22.15 

360 

8 

.27 

25.20 

315 

10 

.12 

29.01 

275 

12 

.04 

32.58 

245 

15 

.01 

37.63 

215 

To  indicate  the  importance  of  considering  the  shape  of  the  life 
density  for  a  product,  operating  characteristic  curves  for  this  plan 
have  been  computed  and  plotted  in  Figure  3  for  other  selected  values  for 
P.  Included  is  a  curve  for  the  case  in  which  p  equals  1.  This  repre¬ 
sents  the  exponential  distribution  widely  used  as  a  model  in  reliability 
and  life -testing  sampling  inspection.  From  these  cuirves  it  may  be  noted 
that  if  the  underlying  distribution  is  actually  of  the  Weibull  form  and 
the  exponential  is  assumed,  the  actual  operating  characteristics  of  the 
plan  may  differ  very  much  from  those  contemplated.  A  discussion  of  the 
sensitivity  of  statistical  procedures  in  current  use  to  departures  from 
the  assumed  exponentiality  will  be  found  in  a  paper  by  Zelen  and  Danne- 
miller,^^ 

It  may  be  noted  in  connection  with  this  example  that  the  MIL-STD-IO^B 
plans  include  matching  double  and  multiple  sampling  plans.  These  offer 
alternative  possibilities  for  reliability  and  life  testing  applications. 

If  incoming  lots  are  either  quite  good  or  quite  bad  (as  is  commonly  the 
case),  substantial  reductions  in  the  number  of  items  that  must  be  tested 
may  be  made.  If  items  are  expensive  and  if  testing  is  destructive  (as 
it  most  likely  will  be  in  life  testing),  a  reduction  in  average  sample 
size  may  be  of  importance.  If  the  test  period,  t,  is  relatively  long, 
however,  the  elapsed  time  required  for  testing  a  second  sample  (or 
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subsequent  ones  in  multiple  sampling)  when  such  samples  are  required  to 
reach  a  decision  may  raise  difficulties. 

Example  (2) 

For  a  second  example  consider  the  case  of  a  manufacturer  who 
knows  that  his  current  production  of  a  certain  component  has  a  mean 
life  of  approximately  52,000  cycles.  Furthermore,  he  has  learned  from 
his  past  experience  with  life  testing  of  these  components  that  he  can 
assume  a  P  value  of  1/2.  A  life  test  period  of  1000  cycles  seems 
justifiable  and  facilities  are  available  for  testing  a  sample  of  150 
Items  from  each  lot.  This  manufacturer  would  like  to  know  what  accept¬ 
ance  criteria  to  apply  so  that  virtually  all  lots  will  be  passed  as 
long  as  the  expected  mean  life  of  52,000  cycles  is  maintained.  He 
would  also  like  to  know  what  consumer  protection  will  be  afforded.  A 
final  question  is  whether  for  this  application  changing  to  a  proposed 
test  period  of  300  cycles  and  a  sample  size  of  $00  items  would  yield 
comparable  or  better  quality  assurance. 

The  first  step  toward  answers  to  these  questions  is  to  compute 
the  (t/p)  X  100  ratio  at  the  mean  life  considered  acceptable.  This 
ratio  is  (1000/52,000)  x  100  or  1.93.  Entering  Table  2  with  this 
value  gives  (with  rough  interpolation)  a  value  for  pi  of  18^  for  a 
0  value  of  1/2.  Assuming  a  probability  of  acceptance  of  .95  is  de¬ 
sired  for  lots  at  the  acceptable  quality  level  of  52,000  cycles  for 
the  lot  mean,  entering  a  table  of  the  cumulative  binomial  distribution 
indicates  an  acceptance  number,  c,  of  3$  items  gives  this  probability 

for  a  sample  size  of  150  items.  This,  then  is  the  desired  acceptance 
criteria. 
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A  simple  measure  of  consumer's  protection  is  to  find  the  lot 
mean  value  at  which  lots  will  like^  be  rejected.  Suppose  a  probability 
of  rejection  of  .90  (of  acceptance  of  .10)  seems  to  be  a  meaningful 
figure.  Reference  again  to  a  binomial  table  indicates  that  for  n  -  150 
and  c  -  35,  the  probability  of  rejection  is  .90  at  a  p«  of  approxi¬ 
mately  26.h%.  Entering  Table  2  with  this  value  gives  a  (t/p)  x  100 
ratio  of  approximately  5.7.  Substituting  a  value  of  1000  cycles  for 
t  in  this  ratio  and  solving  for  p  gives  a  lot  mean  of  17,500  cycles. 

This  figure  for  consumer's  protection  can  be  interpreted  as  follows- 
since  this  quality  -  17,500)  corresponds  to  a  P(A)  -  .10,  under 
this  sampling  plan  (n  -  150,  c  -  35)  on  the  average  90^  of  the  lots 
passed  to  the  consumer  will  have  a  mean  life  of  no  less  than  17,500 
cycles.  This  may  or  may  not  represent  adequate  consumer  protection. 

If  It  does  not,  a  plan  with  a  larger  sample  size  must  be  designed 
and  used. 

An  answer  to  the  third  question  may  be  found  by  making  similar 
computations  for  an  n  of  500  items  and  a  value  for  t  of  300  cycles. 

In  this  case  (t/p)  x  100  equals  (300/52,000)  x  100  or  .56.  From  Table  2 
It  is  found  that  p»  is  approximately  10^  at  this  truncation  ratio. 

Scanning  a  binomial  table  indicates  an  acceptance  number  of  62  will 
give  a  probability  of  acceptance  of  .95  or  more  when  the  sample  size 
is  500  items.  With  this  sample  size  and  acceptance  number,  the  proba¬ 
bility  of  rejection  is  .90  at  a  p'  value  of  approximately  Ihfo,  With 
this  value  for  p',  a  (t/p)  x  100  value  of  approximately  I.15  is  found 
from  Table  1.  Substituting  3OO  cycles  for  t  in  this  ratio  gives  a 
lot  mean  value  of  26,100  cycles  as  compared  to  17,500  cycles  for  the 
first  plan.  Thus  this  combination  of  sample  size  and  length  of  test 
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period  gives  better  discrimination  between  good  and  bad  lots  and  the 
consumer  is  therefore  better  protected. 


cjcample  ( 3  ^ 

In  this  example,  reference  will  be  made  to  a  case  for  which  the 
oemiponent  life  can  best  be  characterized  by  a  mixture  of  two  Wei- 
bul  distributions.  Kao^  gives  an  example  of  this  for  the  life  of 
electron  tubes.  From  the  electron  tube  life  experience,  the  wearout 
failures,  i.e.,  drift  of  electrical  properties  beyond  some  set  limits, 
invariably  occur  near  the  latter  part  of  life.  Hence  the  failures  of 
electron  tubes  are  classified  both  as  of  the  wearout  type  and  as  of  the 
non-wearout  or  catastrophic  type,  each  type  being  represented  by  a  sub- 
population  Of  the  whole.  In  electronic  terms,  these  failure  types  are 
referred  to  as  electrical  rejects  and  inoperative  rejects  respectively. 
The  catastrophic  (or  inoperative  rejects)  sub-population  is  assumed 
to  start  at  time  zero,  i.a,  the  location  parameter  -  o,  when  the 
components  are  exposed  to  risks.  The  wearout  (or  electrical  rejects) 
sub-population  is  assumed  not  to  start  until  some  delayed  period  has 
elapsed,  i.e.,  73  >  0,  since  the  limits  set  on  the  component  drift 
depend  on  many  factors  such  as  environmental  stress,  maintenance  po¬ 
licy,  legal  regulations,  etc.  Since,  in  general,  failures  due  to  wear¬ 
out  and  non-wearout  reasons  are  identifiable,  it  is  possible  to  treat 
the  two  sub-populations  separately. 


Suppose  that  for  some  application  of  electron  tubes,  the  manu¬ 
facturer's  past  experience  indicates  that  the  Weibull  shape  parameter 
associated  with  the  catastrophic  sub-population,  =  1/2  and  that 
associated  with  the  wearout  sub-population  -  3  1/3  are  reasonable 
values,  and  furthermore  that  electrical  drift  or  wearout  failure  has 
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neT,r  been  experienced  prior  to  1000  hoore  of  lif,  (  .  o.  .  loooi 

Suppose  further  that  the  nanufaoturer  knows  that  approxinately  2  V2* 

Of  the  total  tube  faUure.  are  of  the  inoperative  type  and  that  the 
mean  life  of  his  ciirrent  production  is, 

P  -  (.025)  (25,000)  »  (.975)  (11,000)  .  11,325  bra 


-here  .  25,000  is  the  mean  life  of  the  oatastrophio  sub-population 
und  pj  .  11,000  is  that  of  the  wearout  sub-population.  (See  the  ap¬ 
pendix  of  the  paper  by  Kao^  for  the  derivation  of  this  fomwla.) 

A  life  test  period  of  500  hours  for  inoperatives  and  of  5000  hours 
for  elaotrioal  drifts  are  reoomirended  and  acceptance  numbers  Oj^  -  2 
and  (=2  ■  2  for  each  failure  type  are  considered  satisfactory.  What 
are  the  necessary  sample  sloes  so  that  the  producer's  risk  is  no  more 

than  5*2  Also,  what  would  be  the  consumer's  protection  under  this 
sampling  plan? 


To  answer  these  questions,  the  two  sub-populations  are  treated  se 
parately  and  are  denoted  by  subscripts  1  and  2  as  done  before  for  the 
inoperatives  and  electrical  drifts  respectively.  For  inoperatives, 

^  *  AOO)  /25,0OO  ■  2.0,  Entering  Table  2  with 

this  value  gives  a  value  for  p'  of  16.13*  for  a  e  value  of  1/2.  From 
a  binomial  table  with  P(A)  k  .95  and  pj  .  16.13*,  a  value  for  n^  .  5 
is  obtained.  The  same  binomial  table  for  =  5,  =  2  and 

P(A)  a  .10  gives  p.  .  75*.  Entering  Table  2  with  this  value  gi^s 
(tl/dl)  X  100  .  96.5  and  .  (500  x  100)  M.;  .  hours,  a 

value  which  will  be  commented  on  later.  For  electrical  drifts, 

.mist  be  subtracted  from  t^  and  Pj  giving  new  values  for  t^  and  p/ 
equal  to  U,000  and  10,000  rospootively.  Hence,  (tj/p^)  x  100  . 

(11,000  X  100)  /10,000  .  1,0.0  .  Entering  Table  2  with  this  value 
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gives  for  a  figure  of  3.25^  for  a  p  value  of  3  1/3.  Prom  a 
binomial  table  with  P(A)  ^  .95  and  p.  .  3.25^,  it  is  found  that 
ng  -  25,  The  same  binomial  table  for  ng  -  25,  Og  -  2,  P(a)  s  .10 
gi^es  p.  .  20$,  Entering  Table  1  with  this  value  gives  (t/jx^)  x  100  ^ 
71.02,.  Thus  P2  -  (1*,000  X  100)  /71.0U  -  5,631  hours,  which  upon  re- 
adding  72  gives  6,631  hours.  Combining  this  value  of  corrected 
with  obtained  for  inoperatives  gives  the  consumer's  protection^ 
expressed  in  terms  of  a  mean  value  equal  to, 

M.  -  (.025)  (518.21,)  +  (.975)  (6,631)  -  6,2,76  hours. 

This  means  under  the  sampling  plan  of  running  a  life  test  for  inopera- 
tives  Of  n  -  5  and  c  .  2  for  500  hours  and  another  life  test  for 
electrical  drifts  of  n  -  25  and  c  -  2  for  5000  hours,  90$  of  the 
lots  passed  to  consumers  will  have  a  mean  life  of  at  least  6,1,78  hours. 
To  illustrate  the  danger  of  extrapolation  in  a  mixed  distribution 
case,  assume  that  the  second  test  of  5000  hours  duration  was  not  run 
at  all,  then  the  producer  could  only  base  his  conclusion  upon  the 
500-hour  test  and  claim  as  consumers '  protection,  with  90$  confi¬ 
dence,  a  mean  life  of  at  least  .5i8,2l,hour5,  a  result  which  is 
altogether  too  modest. 


REUTIONSHIP  BETl\2EEN  ACCEPTABLE  AMD  UNACCEPTABLE  LOT  QUALITY 

SanipUng  plana  are  most  conveniently  cataloged,  selected,  or  de¬ 
signed  in  toms  of  a  producer's  risk  and  a  consumer's  risk.  Some  lot 
quality  figure  will  be  specified  as  satisfactoiy  and  for  lots  of  this 
quality  or  better  the  probability  of  acceptance  should  be  high,  conven¬ 
tionally  .95  or  more  (the  producer's  risk  of  rejection  small,  .05  or  less) 
For  these  plans  for  Ufe  testing,  this  specification  will  be  a  lot  mean 
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life,  which  P(A)  a  .95.  Likewise,  an  unsatisfactory  quality  level 

will  be  specified  for  which  the  probability  of  acceptance  will  be  low, 
conventionally  .10.  This  specification  will  be  a  lot  mean  life,  n 
at  which  P(A)  s  .10.  (For  other  values  of  producer's  and  consumer's 
risks  reference  may  be  made  to  a  paper  by  Kao.^^) 

In  plan  selection  or  design,  one  objective  is  to  find  a  combination 
of  sample  size  and  acceptance  number  which  simultaneously  yields  the 
desired  values  for  both  the  consumer  risk  and  the  producer  risk.  If 
working  from  tables  of  plans,  the  values  for  lot  quality  at  the  two  risk 
figures  may  be  listed.  In  the  design  of  a  plan,  one  may  cut  and  try 
until  a  suitable  plan  is  found  in  a  manner  suggested  in  Example  2.  Also, 
factors  are  available  which,  in  conjunction  with  the  conversion  tables 
supplied  here,  enable  a  direct  determination  to  be  made.^^ 

A  simple  alternative  solution  for  the  form  of  plan  discussed  here  is 
to  make  use  of  one  of  its  properties,  namely  that  for  a  given  acceptance 
number,  c,  (and  for  a  given  value  forp  )  the  ratio  between  the  lot  means 
at  the  two  risk  values  is  approximately  constant  for  all  values  of  sample 
size,  n.  These  ratios  (or  multipliers)  have  been  determined  for  values 
for  c  ranging  from  0  to  1$  for  each  of  the  various  values  for  3  .  They  are 
presented  in  Table  3.  The  table  values  are  in  the  form  of  multipliers 
for  finding  given  or  using  the  reciprocal  of  the  multiplier, 

for  finding  given  That  is,  (for  which  P(a)  -  .95)  is 

equal  (approximately)  to  p^^^  (for  which  P(A)  =  .10)  times  the  appro¬ 
priate  table  multiplier.  These  multipliers  may  be  used  both  to  assist 
in  evaluating  the  operating  characteristics  of  some  given  plan  and  to 

assist  in  the  design  of  a  plan  to  meet  some  acceptance -inspection  re¬ 
quirement. 
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Example  (lx) 


For  a  certain  purchased  component  the  lot  mean  life  should  be 
at  least  ii,000  hours;  this  value  is  accordingly  chosen  for  p, 

Also,  the  producer  has  been  informed  that  lots  whose  mean  life  is 
10,000  hours  or  more  are  reasonably  sure  of  acceptance  through  the 
sampling  procedure.  Accordingly,  this  value  is  to  be  used  for  p 
A  value  for  p  of  1  can  be  assumed.  A  testing  period,  t,  of  200  hours 
has  been  specified.  Values  for  sample  size,  n,  and  acceptance  number, 
c,  must  be  found  to  meet  these  requirements. 

The  ratio  between  the  two  lot  means,  ,  is  10,000/1*, 000  or 

2=5  .  Examination  of  the  table  of  mean  life  multipliers.  Table  3, 

under  the  column  f or  p  -  1  indicates  that  an  acceptance  number,  c, 

of  10  items  will  give  this  ratio.  The  (t/p)  x  100  ratio  at  p 

is  (200/it, 000)  X  100  or  5.  Entering  Table  2,  the  table  of  p«,  with 

this  truncation  ratio  value  of  5,  gives  a  p'  of  it. 88^.  Reference 

to  a  table  of  the  cumulative  binomial  distribution  or  use  of  the 

Poisson  approximation  for  c  -  10  and  p'  »  .01*88  at  P(A)  -  .10  shows 

that  a  sample  size,  n,  of  315  items  meets  the  requirements.  A  check 

for  this  solution  can  be  made,  if  desired.  For  n  =  315,  c  -  10,  and 

P(A)  =  .95  the  Poisson  approximation  indicates  a  pi  of  1.96^.  Entering 

Table  1,  the  table  for  per  cent  truncation,  with  this  value  for  p', 

a  value  for  (t/p)  x  100  of  approximately  2.0  is  found.  Solving  for 

11.95  yields  (200/p)  x  100  -  2.0  or  p^^^  -  10,000  which  is  the  desired 
value. 

TABLES  OF  SAfIPLING  PLANS 

A  set  of  tables  of  sampling  inspection  plans  has  been  prepared,  one 
table  for  each  of  the  nine  values  for  p  for  which  the  relationship 
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are  presented 


between  p«  and  (t^)  x  100  has  been  established.  These 
as  Tables  i^a  through 

tebl.  ll.t.  «l„es  for  the  aooeptaooe  number,  o,  and  for  the 
numpl.  Slue,  n,  for  a  «rl,tf  of  objeotlva  t/p  ratioa.  The 
Plana  are  daalened  ao  that  if  100  tl»a  the  »tio  between  the  teat  time, 
t,  and  the  mean  life  value  for  the  lot,  p,  „r  (t/p)  x  100  ia  equal  to  or 
greater  than  the  aeleoted  oolumn  value  in  the  table,  the  probabUity  of 
ecoeptanoe,  P(A),  wiu  be  .10  or  leaa.  Stated  otherwiae,  the  plana 
aaaure  with  90?  confidenoe  or  more  the  aooeptanoe  of  lota  for  which  the 
(Vp)  m  100  ratio  ia  equal  to  or  leaa  than  the  aeleoted  oolumn  or  objeo- 
tive  value.  The  ratioa  in  the  oolumn  headinga  (for  which  the  plana  have 

been  deaigned)  may  thua  be  oonaidered  in  the  aame  way  aa  lo*  tolerance 

per  cent  defective  (LTPD)  value*?  oya  ^v,  a 

liairu;  values  are  in  describing  operating  characteris- 

tioa  Of  the  widely  uaed  attribute  «,d  variablea  aooeptanoe  plana. 

It  haa  been  aaaumed  that  in  aooeptanoe  Inapeotion  for  reliability 
the  oonaumer.a  rlak  will  be  of  primary  oonoern.  For  thia  reaaon,  these 
plana  have  been  cataloged  by  P(A)  .  .10  ratioa  which  measure  oonaumer 
protection.  However,  in  addition  for  each  plan  the  (t/p)  x  100  ratio 
ia  given  for  which  the  probability  of  aooeptanoe  ia  .95  or  more.  Each 
such  P(a)  .  .95  ratio  value  ia  shown  in  parentheaea  under  the  corresponding 
aample  site  number.  These  ratio  values  may  be  oonsideted  similar  to 
acceptable  quality  level  (AQL)  values  in  indicating  the  producer's  risk. 

If  the  mean  life  for  the  items  in  the  lot  is  such  that  the  t/p  ratio  is 
equal  to  or  leas  than  the  tabulated  value,  there  is  assurance  with  con¬ 
fidence  of  95*  or  more  that  the  lot  will  be  accepted. 

The  two  ratio  values,  one  in  the  oolumn  heading  and  the  other  in 
parentheaea  below  the  aample  aiae  number,  broadly  describe  the  operating 
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characteristics  of  each  plan  and  so  form  a  basis  for  making  an  appro¬ 
priate  choice  for  any  acceptance  inspection  application.  Those  values 
may  also  be  used  to  determine  approximately  the  operating  characteristics 
of  any  acceptance  plan  that  has  been  specified  or  that  is  in  use  and 
for  which  n  and  c  match  closely  one  of  the  table  plans.  It  is  easy  to 
convert  these  ratios  to  hours,  cycles  or  some  other  measure  of  lifelength 
to  fit  the  product  and  test  specifications  involved.  This  will  be  illus¬ 
trated  by  two  examples  which  follow  later. 

In  the  preparation  of  these  plans,  binomial  tables  prepared  by 
17 

Grubbs  were  employed  for  values  for  c  up  to  9  and  for  n  up  to  1^0. 

For  higher  values  of  c  and  for  values  for  n  up  to  60  or  so,  the  Pearson 
tables  of  the  incomplete  beta-function  were  used.^®  Higher  values  of  n 
were  deterinined  by  the  Poisson  approximation,  using  a  table  of  npi  ‘ 
values  prepared  by  Cameron.^^  The  Poisson  match  was  checked  and  was 
found  close,  even  for  the  smaller  sample  sizes  and  large  values  for  p*. 
The  slight  differences  that  may  exist  in  some  cases  is  on  the  conserva¬ 
tive  sidej  the  value  for  n  is  slightly  larger  than  that  theoretically 
required.  As  this  is  primarily  an  exploratory  study,  plans  showing 
extremely  large  sample  sizes  have  been  included  to  indicate  the  order 
of  magnitude  involved  and  not  with  the  expectation  that  samples  of  this 
size  would  ordinarily  be  used. 

Example  (5) 

An  acceptance  inspection  plan  is  required  which  will  assure 
with  90%  confidence  a  mean  life  for  items  of  ij,000  hours  or  more  for 
each  lot  accepted.  Also,  it  will  be  desirable  to  assure  the  pro¬ 
ducer  that  if  the  mean  life  for  items  in  a  lot  is  25,000  hours  or 
more,  there  will  be  a  high  probability  (.95)  of  its  acceptance. 
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A  test  period  of  I4OO  hours  for  the  inspection  of  sample  items  has 
been  specified.  Through  past  experience  it  has  been  determined  that 
the  distribution  of  item  life  is  of  the  Weibull  form  with  3  equal 
to  approximately  1/2. 

For  these  sampling  plan  specifications  100  times  the  ratio  of  test 
time,  t,  to  mean  life,  is  (UOOA^OOO)  x  100  or  10  for  which  a 

probability  of  acceptance  of  .10  or  less  is  desired.  At  the  .9? 
probability  value  the  ratio  is  (iiOO/25,000)  x  100  or  1.6.  A  plan 
approximating  this  may  be  found  in  Table  J|b  which  gives  plans  for 
distributions  for  which  P  -  1/2.  The  column  for  which  (t/p)  x  100  -  10 
is  entered  and  scanned  for  the  ratio  value  1.6  among  the  values  listed 
in  parentheses.  This  value  is  found  well  down  in  the  column.  The 
corresponding  sampling  inspection  plan  specifies  a  sample  size,  n,  of 
h3  and  an  acceptance  number,  c,  of  11. 

Excunple  (6^ 

A  sampling  inspection  plan  specifies  that  a  random  sample  of 
3000  items  be  drawn  from  the  lot  and  tested  for  a  period  of  h80 
hours.  If  no  more  than  7  items  fail  before  the  end  of  the  test 
period,  the  lot  is  to  be  accepted;  if  more  than  7  items  do  not  live 
through  the  test  period,  the  lot  is  to  be  rejected.  Life  measure¬ 
ments  for  past  inspection  and  research  for  the  product  to  which  the 
plan  is  to  be  applied  indicate  the  distribution  is  of  the  Weibull  form 
with  p  equal  to  approximately  1  2/3.  The  prospective  user  of  this  plan 
would  like  to  know  what  quality  protection  will  be  given.  Inspection 
of  Table  hd  which  lists  plans  for  p  -  1  2/3  discloses  a  plan  match¬ 
ing  reasonably  well  the  one  specified,  the  plan  for  which  c,  the 
acceptance  number,  is  7  and  n,  the  sample  size,  is  3,019.  For  this 


table  plan  the  (t/p)  x  100  ratio  at  P(A)  -  .10  is  i*.  Substitution  of 
the  specified  test  period  length  of  2|80  hours  for  t  gives  (U80/p)  x 
100  ■  ii.  Solving  for  p  gives  12,000  hours  as  the  mean  value  for  item 
life  for  the  lot  for  which  the  probability  of  acceptance  is  .10  or 
less.  A  similar  substitution  for  t  using  the  ratio  at  which  P(A)  ■ 

.95  gives  (h60/ii)  x  100  -  2.1.  Solving  for  p  again  gives  23,000 
hours  as  a  lot  mean  value  for  which  the  probability  of  acceptance  is 
.95  or  better.  The  values  for  the  lot  mean  at  these  two  probability 
values  broadly,  but  very  practically,  describe  the  operating  charac¬ 
teristics  of  the  specified  plan. 

In  the  use  of  these  tables  of  plans,  several  points  of  practice 
should  be  noted.  First,  in  using  the  p»  values  associated  with  values 
for  (t/p)  xlOO  for  the  V/eibull  distribution  to  find  values  for  c  and  n, 
the  binomial  probability  distribution  has  been  used.  Hence  the  size  of  the 
lot  should  be  relatively  large  compared  to  the  size  of  the  sample  for 
the  stated  probability  values  to  precisely  apply.  Second,  if  a  plan  is  not 
available  for  which  a  (t/p)  x  100  ratio  in  the  column  headings  matches 
closely  the  desired  ratio,  to  be  conservative,  a  plan  should  be  chosen 
from  the  column  with  the  next  smaller  ratio  heading.  This  will  assure 
with  confidence  greater  than  90^  the  specified  mean  life  for  acceptance. 

If  the  acceptable  quality  level  (the  ratio  or  mean  life  for  which  P-(A)=.95) 
must  also  be  guaranteed  and  a  matching  ratio  value  is  not  found  in  the 
selected  column  of  plans,  a  plan  with  the  next  greater  value  should  be 
selected.  Lots  equal  to  or  better  than  the  specified  "acceptable  quality" 
will  have  an  assurance  of  greater  than  95^  of  being  accepted.  With 
proper  care,  some  rough  interpolation  may  be  employed  between  listed 
sample  sizes  (either  down  or  across  the  table  or  both)  to  find  a  new 
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plan  having  more  nearly  the  desired  characteristics.  Finally,  if  a 
plan  is  found  for  which  the  desired  and  given  ratios  closely  match  but 
for  practical  reasons  it  seems  desirable  to  round  off  the  sample  size  to 
the  nearest  number  ending  in  zero  or  five,  such  rounding  off  should  be 
done  to  the  next  larger  size.  This  will  assure  retention  of  the  proba¬ 
bility  values  of  .10  or  less  for  the  ratios  given  in  the  column  headings. 
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Figure  I.  Plot  of  the  Wei  bull  Probability  Density  Function  for  Various  Values  of  P  -  (t=I) 
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TABLE  1 


Table  of  Values  for  Per  cent  Truncation,  (t/jx)  x  100 


.76 

2.22 

5.55 

19.28 

69.48 


.030 

.040 

.050 

.065 

.080 

.86 

1.02 

1.18 

1.37 

1.56 

1.95 

2.26 

2.53 

2.88 

3.19 

4.40 

4.93 

5.39 

5.98 

6.50 

9.77 

10.65 

11.40 

12.32 

13.13 

14.52 

15.60 

16.49 

17.62 

18.56 

21.50 

22.77 

23.82 

25.10 

26.16 

.10 

.12 

.15 

.20 

.25 

1.78 

1.98 

2.26 

2.69 

3.08 

3.57 

3.92 

4.37 

5.07 

5.64 

7.U 

7.65 

8.36 

9.39 

10.27 

14.03 

14.82 

15.84 

17.27 

18.47 

19.62 

20.53 

21,71 

23.33 

24.68 

27.36 

28.37 
29.67 
31.43 
32.87 

.001 

.002 

.003 

.30 

.40 

.50 

.65 

.80 

3.44 

4.07 

4.67 

5.46 

6.19 

6.18 

7.14 

7.99 

9.12 

10.11 

11.05 

12.39 

13.55 

15.06 

16.36 

19.51 

21.27 

22.75 

24.62 

26.21 

25.83 

27.76 

29.36 

31.35 

33.03 

34.09 

36.12 

37.76 

39.81 

41.50 

.005 

.007 

.011 

.020 

.032 

1.01 

1.21 

1.51 

2.02 

2.53 

7.08 

7.90 

9.07 

10.77 

12.33 

11.31 

12.40 

13.87 

16.03 

17.95 

17.90 

19.26 

21.08 

23.67 

25.90 

28.03 

29.62 

31.68 

34.56 

36.98 

34.93 

36.57 

38.68 

41.59 

44.01 

43.40 

45.02 

47.09 

49.90 

52.21 

.047 

.083 

.13 

.23 

.35 

3.05 

4.08 

5.13 

6.72 

8.34 

13.78 

16.42 

18.84 

22.15 

25.20 

19.69 

22.79 

25.58 
29.25 

32.59 

27.89 

31.35 

34.35 

38.28 

41.72 

39.09 

42.69 

45.71 

49.57 

52.88 

46.09 

49.59 

52.50 

56.18 

59.29 

54.17 

57.45 

60.13 

63.46 
66.26 

.56 

.82 

1.32 

2.49 

4.14 

10.54 

12.78 

16.25 

22.31 

28.77 

29.01 

32.58 

37.63 

45.51 

52.99 

36.63 

40.34 

45.48 

53.30 

60.53 

45.82 

49.50 

54.49 

61.85 

68.47 

56.73 

60.11 

64.60 

71.04 

76.67 

62.85 

65.96 

70.05 

75.83 

80.80 

69.44 

72.18 

75.73 

80.68 

84.89 

35.37 

51.08 

69.31 

104.98 

160.9)1 


60.29 

74.79 

89.32 

115.23 


67.39 

80.64 

93.95 

115.61 

il.'a  T  I. 


74.62  81.79 

86.15  91.09 

97.33  99.82 

114.92  113.06 


T  oz:  I 


85.26  88.62 

93.27  95.22 

100.67  101.21 

111.68  109.98 


TABLE  2 


Table  or  Probability  Valaaa  .t  Tnmoation  Polat,  p.  (« 


(t/ji)x  100 


Shape  Parameter  -  p 


1/3 

.010 

8.09 

.012  1 

8.57 

.015 

9.20 

.020 

10.08 

.025 

10.82 

.030 

11.1*5 

.01*0 

12.53 

.050 

13.1*3 

.065 

11*.56 

.080 

15.52 

.100 

16.61 

17.56  U.78 
13*78  5.33 
20. h6  6.13 
21.86  6.83 

23.06  7.h5 
25.06  8.56 
26.71  9.52 
28,76  10.78 
30.U7  11.88 

32. ho  13.19 
3h.03  lh.3B 
36.12  15.90 
38.91;  18.13 
hi. 22  20. Oh 


1/2 

1 

1.1*0 

.010 

1.51* 

.012 

1.72 

.015 

1.98 

.020 

2.21 

.025 

2.1*2 

.030 

2.79 

.01*0 

3.11 

.050 

3.51* 

.065 

3.92 

.080 

h.37 

.10 

h.78 

.12 

5.33 

.15 

6.13 

.20 

6.83 

.25 

1,00 

1.19 

1.1*9 

1.98 

2.1*7 


1*3.11* 

1*6.28 


21.73 

21*.61* 


1*8.80  27. U 
50.90  30.27 
51*.  30  32,97 


1  2/3 


2  1/2  3  3/3 


10.0 

12 

15 

20 

25 


56.98 

59.19 

61.92 

65.1*5 


36.06 

38.73 

1*2.17 

1*6.87 


68.17  50.69 

70.37  53.91 
73.79  59.12 


76.36 

79.28 

81.1*9 

83.75 


63.21 

68.02 

71.77 

75.69 


9.52 

11.31 

13.93 

18.13 

22.12 

25.92 

32.97 

39.35 

1*7.80 

55.07 

63.21 


.052 

.076 

.12 

.18 

.011 

.018 

.031 

.01*9 

.001 

.002 

.001* 

.007 

.21* 

.39 

.56 

.89 

1.22 

.071 

.13 

.20 

.33 

.50 

.012 

.021* 

.Ola 

.080 

.13 

.001 

.002 

.003 

.008 

.015 

.001 

.003 

1.77 

2.39 

3.1*5 

5.51 

7.89 

.78 

1.12 

1.75 

3.09 

1*.79 

.23 

.37 

.61* 

1.32 

2.29 

.033 

.060 

.13 

.33 

.69 

.007 

.011* 

.031* 

.11 

.26 

10.56 

16.1*7 

22.98 

33.26 

1*3.51* 

56.35 

6.82 
11.81 
17.83 
28.21* 
39.51 
51*. 1*1 

3.59 

7.23 

11.28 

22.32 

31*.59 

52.36 

1.26 
3.25 
6.72 
15.37 
28.35 
50.  la 

.55 

1.71 

1*.13 

11.35 

21*.15 

1*9.08 

.16 

.67 

2.02 

7.29 

19.25 

1*7.93 


TABLE  ka. 

Table  of  San9>ling  Plans  for  P  ■  3/3 


(t/ti)  X  100  Ratio  for  which  P(A)  -  .10  (or  Less) 
^ _ 


:  I  1 - 1 


(•ou)  (.02) 


(t/n)  X  100  ratios  In  parentheses  are  for  P(A)  •  .95  (or  more) 


TABLE  hb 

Table  of  Sampling  Plans  for  P  -  1/2 


(t/n)  X  100  ratios  in  parentheses  are  for  P(A)  ■  .95  (or  more) 


TABLE  Uc 

Table  of  Sampling  Plans  for  P  •  1 


(t/p)  X  100  Ratio  for  which  P(a)  -  .10  (or  less) 


2.5 

1 

92 

231 

(.06) 

(.02) 

158 

389 

(.22) 

(.09) 

216 

533 

(.38) 

(.15) 

271 

669 

(.50) 

(.20) 

•5  0.25  0,1  0.05  0.025  0.01 


9,212  230-2 
156-2  389-2 
213-2  532-2 
267-2  668-2 
800-2 


109 

(3.0) 

216 

(1.5) 

1|27 

(.77) 

l,051i 

(.31) 

122 

(3.3) 

2U2 

(1.7) 

k77 

(.83) 

1,178 

(.3li) 

135 

(3.5) 

267 

(1.8) 

527 

(.89) 

1,300 

(.36) 

1U7 

(3.7) 

292 

(1.9) 

576 

(.9W 

l,ii21 

(.38) 

162 

(3.9) 

316 

(2.0) 

62U 

(1.0) 

1,5U1 

(.liO) 

211-2 

li21-2 

105-3 

(.02) 

(.01) 

235-2 

li71-2 

U8-3 

(.02) 

(.01) 

260-2 

520-2 

130-3 

(.02) 

(.01) 

281i-2 

568-2 

1U2-3 

(.02) 

(.01) 

308-2 

616-2 

15U-3 

(.02) 

(.01) 

3,320 

6 ,640 

166-2 

(.21) 

(.10) 

(.04) 

3,557 

7,113 

178-2 

(.22) 

(.11) 

(.04) 

3,792 

7,584 

190-2 

(.22) 

(.11) 

(.05) 

4,026 

8,052 

'  201-2 

(.23) 

(.12) 

(.05) 

4,260 

8,517 

213-2 

(.24) 

(.12) 

(.05) 

(t/n)  X  100  ratios  in  parentheses  are  for  P(A)  -  .95  (or  more) 


^  j  j  ^^8*1**  following  the  dash  in  sample  size  numbers  shows  the  number  of  zeros  to 
addj  for  exanq>le,  l5k-3  •  l5U,000. 


TABLE  hd 

Table  of  San^^Ung  Plans  for  P  ■  1  2/3 


(t/p.)  X  100  ratios  in  parentheses  are  for  P(A)  ■  ,95  (or  more) 

The  figure  following  the  dash  in  sample  size  numbers  shows  the  number  of 
to  add;  for  example,  218-2  ■  21,800, 


zeros 


IA61E 

Table  of  Saii5>llng  Plana  for  P  -  2 


TABLE  Uf 

Table  of  Sampling  Plans  for  p  «  2  1/2 


(t/n)  X  100  Ratio  for  which  P(A)  ■  .ip  (or 

' '  _  ' 


100 1  50  iio|  2S|  15|  12|  10 


8  6.5 


0  k 
(20) 

1  6 
i(37) 

2  9 

!(h5) 
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The  figure  following  the  dash  in  sample  size  numbers  shows  the 
zeros  to  add;  for  example,  30li-3  -  30^,000. 
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The  figure  following  the  dash  in  sample  size  numbers  shows  the  number  of 
zeros  to  add;  for  exanqple,  203-2  «  20,300. 
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